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ABSTRACT 

Distance measurements to Type la supernovae (SNe la) indicate that the Universe is 
accelerating and that two-thirds of the critical energy density exists in a dark-energy com¬ 
ponent with negative pressure. Distance measurements to SNe la can be used to distinguish 
between different possibilities for the dark energy, and if it is an evolving scalar held, to 
reconstruct the scalar-held potential. We derive the reconstruction equations and address 
the feasibility of this approach by Monte-Carlo simulation. 


1 Introduction 


There is now prima facie evidence that the Universe is flat and that the critical energy density 
is 1/3 matter and 2/3 something else with large, negative pressnre. The simplest possibility 
for the latter component is vacnnm energy (cosmological constant) |l|; other possibilities 
include a frustrated network of topological defects 0 and an evolving scalar field 0, ||, 
called quintessence by the authors of Ref. |^. All have effective bulk pressure that is very 
negative, p < —p/3; for the cosmological constant p = —p and for a frustrated defect network 
p = —^p where N is the dimension of the defect. In this paper we discuss the use of SNe la 
to distinguish between these possibilities and to probe the scalar-field potential associated 
with the quintessence field. 

Backing up for a moment, the evidence for flatness comes from measurements of the 
multipole power spectrum of the cosmic background radiation (CBR) which show a peak 
around I ~ 200 as expected for a fiat Universe |^. A variety of dynamical measurements 
of the mean matter density indicate that Qm = 0.4 ±0.1 0]. Recent measurements of 
the distances to more than 50 SNe la out to redshift z ^ 1 indicate that the expansion is 
accelerating rather slowing down If correct, this implies the existence of an unknown 
component to the energy density with pressure px = wxPx ^ —Px/3 that contributes fix ~ 
0.6 1^. This fits neatly with the determinations that Qm ~ 0.4 and Hq (= ± ^m) ~ 1- 

While this accounting is not yet definitive - and could possibly change dramatically - it is 
worth thinking about how to distinguish between the different possibilities suggested for the 
unknown energy component . 

The key difference between quintessence and the other two possibilities is that the effective 
equation of state, wx = Px/px, can vary with time and can take on any value. The 
combination of SNe la measurements and high-precision measurements of the multipole 
power spectrum expected from the MAP and Planck Surveyor satellites may be able to 
discriminate between constant and varying wx [0- If wx is found to vary and/or is not 


equal to —N/3 {N = 1,2,3), the next question is how best to probe the “quintessence 
sector.” While anisotropy of the CBR will be very powerful in determining many important 
cosmological parameters, as we now explain, it has less potential to probe the scalar-field 
potential than SNe la measurements. The fundamental reason is simple: CBR anisotropy 
primarily probes the Universe at redshift z ~ 1000 when the ratio of dark-energy density to 
matter density was tiny (-C 10“®); the SNe la probe the Universe at recent epochs when the 
dark-energy density is beginning to dominate the matter density. 

Quintessence has three basic effects on CBR anisotropy. The most significant is in de¬ 
termining the distance to the last-scattering surface (Robertson-Walker coordinate distance 
to redshift 2 ; ~ 1100), which sets the geometric relationship between angle subtended and 
length scale. However, all models with the same distance to the last-scattering surface will 
have essentially the same multipole power spectrum. The second and third effects break this 
degeneracy, but are less significant and/or powerful; late-time ISW and slight clumping of 
the scalar field (spatial inhomogeneity induced by the lumpiness in the Universe) only affect 
the lower-order multipoles, which can be less well determined because of cosmic variance 
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Supernovae on the other hand may be able to unravel the essence of quintessence. This 
is because accurate supernovae distance measurements can map out r{z) to redshift 2 : ~ 1 or 
perhaps higher, and this is when quintessence is becoming dynamically important and where 
most of the “scalar-held action” is occurring. [The quantity we focus on, coordinate distance 
to redshift z, r^z), is simply related to the quantity measured by observers, luminosity 
distance, di = {I + z)r{z).] Shortly, we will show the fact the scalar-held action occurs at 
modest redshifts is a natural consequence of quintessence. 

In the next Section we will derive the reconstruction equations for the scalar-held po¬ 
tential, and in the following Section we will address the practicality of this approach with 
simulated data and Monte-Carlo realization of reconstruction. We hnish with a brief sum¬ 
mary and concluding remarks. 


2 Reconstruction Equations 


We assume a hat Universe with two components to the energy density: nonrelativistic mat¬ 
ter, which presently contributes fraction VLm to the critical density, and a single, homoge¬ 
neous scalar held 0 (for the problem at hand, its slight clumping can be neglected). The 
fundamental equations governing our cosmological model are: 
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where r{z) is the Robert son-Walker coordinate distance to an object at redshift ; 2 , the matter 
density pm = ^mPctu = + z)^/87iG, prime denotes derivative with respect to 0, 

and the energy density and pressure of the evolving scalar held are: 

P<P = ^ 0 ^ + '' 4 ( 0 ) ( 5 ) 

P<p = ^ 0 ^-'' 4 ( 0 ) ( 6 ) 

Note too: dz/dt = —(1 -f z)H{z) = —(1 -|- z)/{dr/dz) and = p^/pcut = 1 — ^m- Since 
the relative fractions of critical density in matter and quintessence evolve with time, it is 
important to remember that Qm and refer to the present epoch. 
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As an aside, and before deriving the reconstruction equations, we will show why quintessence 
models are likely to predict interesting scalar-held dynamics recently. As an indicator of “in¬ 
teresting” scalar-held dynamics, consider the time derivative of the ehective equation of 
state, ^ 'Pcj)/P(i>i 

d\nw^ ^d\nPE/dhia KE-PE 

dlna “ 1 - KE/PE ~ KE‘^ - PE‘^ ^ ^ 

where KE = PE = V{(j)), and the equality follows from using the equation of motion 
for (j). If quintessence is to be distinguishable from a cosmological constant, then Wfj, must 
diher from —1, which implies that the kinetic and potential terms are comparable. Further, 
barring accidental (or pre-arranged cancellations), Eq. (|^) then implies that dlnw/dlna is 
presently of order unity. 

On to reconstruction. Since r{z) is determined by H{z) and H{z) is a function of the 
scalar held, one should be able to write down equations for V{(j)) and (p in terms of r{z). The 
following is a parametric solution for V{(j)) and dcp/dz, in terms of r[z), dr/dz and d?r/dz‘^\ 
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where the upper (lower) sign applies if 0 > 0 (< 0). The sign in fact is arbitrary, as it can 
be changed by the held redehnition, 0 —0. The actual value of 0 cannot be determined 

by reconstruction: 0 can be shifted by an arbitrary constant with no cosmological ehect (the 
form of the potential of course changes). 

In integrating the reconstruction equations it is useful to dehne dimensionless quantities. 


X = Hat 

f = H^r 

0 = 0/mpi 

H = H / Hq = ^ ffM(l + z)^ + uj + 

a;(0) = V{(j)m-p\) / (“iHl/STTG) 


The diherential equations governing r, 0 and a become 
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and the construction equations are 
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The boundary conditions, expressed at the present epoch, are: ifo = 1, 0 < a;(0o) < 1 — 
and d^/dx = \J^[l - -o;]. 

Finally, without recourse to a scalar-held model for the unknown, negative pressure 
component, one can derive a reconstruction equation for the bulk equation of state, wx = 
Px!Pxi as a function of redshift. The equation of motion for the X-component, 

d\npx = -3(l + n;x)dln(l + ^) ^ px{z) = g^/; (i+-x)din(i+.) ^ 

ottG 


is used in place of Eq. (1), and the reconstruction equation for wx{z) is derived as before. 
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Using Eq. ( 0 ). SNe la measurements alone can be used to determine wx and address its 
time variation. 


3 Simulating Reconstruction 

Here we investigate the feasibility of our approach and estimate the inherent errors by means 
of simulated data and Monte-Carlo realization. Our procedure is straightforward: 

• Pick a potential U(0), matter density Om, and values for 0o and 0o (consistent with 

1 ) 

• Compute the evolution of 0, a{t) and r{z) for this quintessence model by evolving 0(t) 
and a{t) back in time 

• Realize the model by simulating SNe la measurements: for Zi = {i — l)Az and i = 1 
to N, ri = r{zi) + dri where dri is drawn from a Gaussian distribution with zero mean 
and variance uri{z) {a is the relative error in the luminosity distance) 

• Fit the simulated data with a (low-order) polynomial and numerically compute U(0) 
from the reconstruction equations 

• Repeat one thousand times to estimate the error in reconstructing V (0) 

• This procedure actually reconstructs a;(0); to get V{cj)) one simply takes a value for 
Hq] for our results we took Hq = 70 km sec“^ Mpc“^ 
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Figure 1: The 95% confidence interval for the reconstructed potential assuming luminosity 
distance errors of 5% and 2% (shaded areas) and the original potential (heavy line). For this 
reconstruction, VLm = 0.3, %(0) = Vo[l + cos(0//)], Vq = (4.65 x 10“^eV)^, f/mpi = 0.154, 
N = 40, and = 1.0. The simulated data were fit by a fourth-order polynomial in 

Before presenting some results, we should elaborate on a few technical details. In fitting 
a polynomial to r* we have tried third-, fourth-, fifth-, and sixth-order polynomials; all give 
similar results. Because we are taking derivatives of r{z), the use of higher-order polynomials 
only introduces numerical “noise” and is not useful. We have varied the number of redshift 
bins N from 20 to 100, and Zmax = {N — l)Az from 1 to 1.5; the errors in the reconstruction 
scale roughly as expected, 1/'/N. We also tried using a Gaussian distribution in redshift 
(with more data points at small 2 ; than at high z)] the results change little relative to the 
case of linear distribution. 

We have reconstructed several potentials; here, we present results for the exponential 
and cosine potentials considered previously ^ |^. Shown in Figs. 1 and 2 are the original 
potential and the 95% confidence intervals for the reconstructed potential for data with 
5% and 2% luminosity distance errors. The confidence intervals are obtained by requiring 
that 950 of the 1000 Monte Carlo realizations give a value of the potential in that interval. 
The error in the reconstructed potential is mostly attributable to the (Pr/dz^ terms in the 
reconstruction equations, reflecting the fact that it is extremely difficult to infer the second 
derivative of the noisy data. In particular, an uncertainty of 10% in VLm does not change the 
reconstruction confidence regions appreciably. 
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Figure 2: The 95% confidence interval for the reconstructed potential assuming luminosity 
distance errors of 5% and 2% (shaded areas) and the original potential (heavy line). For 
this reconstruction, VLm = 0.4, K(0) = Vq exp(—/30/mpi), Vq = (2.43 x 10“^ eV)"^, P = 8, 
N = 40, and Zmax = 1-5. The simulated data were fit by a fourth-order polynomial in 

Finally, Fig. ^ shows the reconstruction of the equation of state for an exponential po¬ 
tential. Even though the uncertainties in the reconstruction are large, we are still able to 
distinguish this quintessence model from a constant equation of state. In particular, using 
the statistic for r{z) and the simulated data with 2% errors, all constant equations of 
state can be ruled out with 99.9% confidence, and the most plausible constant equations 
of state, wx = —N/3 for N = 0, 1,2,3, can be ruled out with more than 99.9% confi¬ 
dence. (Of course, the ability to discriminate between constant wx and varying wx depends 
upon how Wx varies.) Fig. also nicely illustrates a general feature of reconstruction: Be¬ 
cause the fractional contribution of dark energy rapidly decreases with increasing redshift, 
Px/pm oc (1 -|- , reconstruction beyond redshift of around 2 ; ~ 0.8 becomes extremely 

difficult. Thus, for the purposes of reconstruction, data at redshifts greater than unity are 
of very limited value. 

4 Discussion 

If correct, the discovery that the expansion of the Universe is speeding up rather than 
slowing down is one of the most important discoveries of this century. It implies that the 
primary component of the cosmos today is dark energy with large negative pressure and of 
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Figure 3: The 95% confidence interval for the reconstructed equation of state of the unknown 
component assuming luminosity distance errors of 5% and 2% (shaded areas) and the original 
equation of state (heavy line). For this reconstruction, VLm = 0.4, %(0) = Vq exp(—/90/mpi), 
Vo = (2.43 X 10“^ eV)"^, /? = 15, iV = 40, and = 1-0. The simulated data were fit by a 
third-order polynomial in 2 ;. 


unknown composition. The implications for fundamental physics are equally profound as all 
the possibilities for the dark-energy component are deeply rooted in fundamental physics: 
vacuum energy, a frustrated network of topological defects and an evolving scalar field. 

The measurement of distances to supernovae out to redshift of order unity, which led to 
this discovery, may also be of great utility in determining the character of the dark energy. In 
this paper we have derived the equations that relate the equation of state of the dark-energy 
component and, in the case of quintessence, the scalar-field potential to measurements of 
the luminosity distance. By use of Monte-Carlo simulation we have shown how this might 
be done in practice and addressed the feasibility of this technique. 

While ours is a preliminary investigation, the results indicate that using SN la measure¬ 
ments to probe the dark-energy component is promising. However, important questions and 
issues remain before one can be confident that this technique can be used in practice. Some 
involve the technical details of how our method might be implemented: What is the optimal 
distribution of supernova redshifts for probing wx and H(0)? Would a likelihood analysis 
for a parameterized fit to the potential be more powerful than our reconstruction approach? 
We have already begun to address these questions with some success. For example, we have 
found that Fade approximants are a much better way to represent the observational data 
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than polynomials (splines may do even better); for many potentials a linear distribution in 
redshift minimizes the area of the 95% conhdence region for the reconstructed potential. 

Foremost among the open issues is the reliability of SNe la as distance indicators. Cur¬ 
rently, the distances errors are estimated to be of order 10% to 20% (per supernova) |I2| . 
They arise from a variety of sources: reddening due to the host galaxy or intergalactic dust; 
intrinsic dispersion in the brightness - decline relationship (Phillips relation |^) used to 
calibrate the SNe la; possible systematic (with redshift) evolutionary effects; and possible 
dependence upon the differing chemical composition of the supernova progenitors. There is 
much activity, both theoretical and observational, directed at better understanding type la 
supernovae and so we can hope for improvement in reducing and/or better understanding 
systematic errors. Further, if the intrinsic scatter (both statistical and systematic) in the 
brightnesses of SNe la is random (not correlated with redshift), then distance error can be 
beaten down by more measurements. For example, allowing 10% distance error per super¬ 
nova, 3% errors in r{z) could be obtained if 10 supernovae are measured at each redshift, 
increasing the total needed for our reconstruction method to 500 or so. Finally, the two 
groups have proven that discovering 1000s of SNela over the next decade is a very realistic 
goal, and with a better understanding of SNe la, one could cull a large sample to form a 
smaller, higher quality sample with smaller systematic error and statistical scatter. 

To summarize, distance measurements to type la supernovae with redshifts less than 
order unity have the potential to shed light upon the nature of the dark energy. Based upon 
existing SNe la measurements and their uncertainties, it appears that obtaining data of the 
quantity and quality required to probe the dark-energy component using the method we 
have described, while not guaranteed, is also not unrealistic. 
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